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The influence of the presence of hyperons in dense hadronic matter on the quantum nucleation 
of quark matter is examined at low temperatures relevant to neutron star cores. We calculate the 
equation of state and the composition of matter before and after deconfinement by using a relativistic 
mean-field theory and an MIT bag model, respectively; the case in which hyperons are present in 
the hadronic system is considered, together with the case of the system without hyperons. We And 
that strangeness contained in hyperons acts to reduce a density jump at deconflnement as well as 
a lepton fraction in the hadronic phase. As a result of these reductions, a quark matter droplet 
being in a virtual or real state has its effective mass lightened and its electric charge diminished into 
nearly zero. The Coulomb screening of leptons on the droplet charge, which has significance to the 
droplet growth after nucleation in the absence of hyperons, is thus shown to be of little consequence. 

If the effective droplet mass is small enough to become comparable to the height of the potential 
barrier, the effect of relativity brings about an exponential increase in the rate of droplet formation 
via quantum tunneling, whereas the role played by energy dissipation in decelerating the droplet 
formation, dominant for matter without hyperons, becomes of less importance. Independently of 
the presence of hyperons, the dynamical compressibility of the hadronic phase is unlikely to affect 
the quantum nucleation of quark matter at temperatures found in neutron star interiors. For matter 
with and without hyperons, we estimate the overpressure needed to form the first droplet in the 
star during the compression due to stellar spin-down or mass accretion from a companion star. 

The temperature at which a crossover from the quantum nucleation to the Arrhenius-type thermal 
nucleation takes place is shown to be large compared with the temperature of matter in the core. 

We also determine the range of the bag-model parameters such as the bag constant, the QCD fine 
structure constant, and the strange quark mass where quark matter is expected to occur in the star. 

PACS number(s): 26.60.-l-c, 12.38.Mh, 64.60.Qb, 97.60.Jd 


I. INTRODUCTION 

A neutron star core, at densities near and just above the normal nuclear density no « 0.16 fm“^, consists of uniform 
nuclear matter being electrically neutral and roughly in /3 equilibrium between its components such as neutrons, 
protons, electrons, and muons. For still higher densities, it is possible that various forms of matter including hyperons, 
meson (pion or kaon) condensates, quark matter droplets, and/or quark-gluon plasmas might become energetically 
favorable; their presence is predicted to play an important role in determining the structure and evolution of neutron 
stars [1]. Basically, these forms are thought to occur via pressure-induced phase transitions, whose nature, both static 
and dynamical, is still uncertain in the absence of adequate information about the properties of matter in the core 
such as the mutual interactions between its possible constituents and the equation of state. 

The possible existence of quark matter in neutron stars started to be considered about two decades ago mainly by 
comparing the energies of hadronic and quark matter at zero temperature [2]. Using a separate physical description 
of the two phases, the transition from uniform hadronic matter to uniform quark matter, which proceeds suddenly at 
constant pressure, was expected to occur in the density range ~ S-IOuq. The presence of uniform quark matter would 
enhance neutrino luminosity and neutron star cooling via quark Urea processes [3] as well as soften the equation of 
state effectively [4]. However, the transition in the /3-equilibrated matter, which contains two conserved charges, i.e., 
electric charge and baryon number, should proceed through a mixed phase persisting over a finite range of pressure, 
as predicted by Glendenning [5] according to Gibbs’ criteria for phase equilibrium. In the low pressure regime of the 
mixed phase, as asserted by Glendenning [5] and Heiselberg et al. [6], quark matter droplets form a Goulomb lattice 
embedded in a sea of hadrons and in a roughly uniform sea of electrons and muons. This spatial structure can arise 
because the presence of strange and down quarks in the negatively charged quark phase plays a role in decreasing the 
electron and muon Fermi energies and in increasing the proton fraction in the positively charged hadronic phase. With 
increasing pressure, the shape of quark matter changes from spheres to rods and then to plates, until the part of quark 
matter and that of hadronic matter begin to be replaced by each other. After further changes of the shape of hadronic 
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matter from plates to rods and then to spheres, the system turns into uniform quark matter at the highest pressure. 
Such structure of one- to three-dimensional periodicity, which has a characteristic length ~ 10 fm, would contribute 
not only to quake phenomena, possibly relevant to pulsar glitches, via its elastic properties [7], but also to enhanced 
neutron star cooling due to opening of the quark phase space for neutrino-generating processes that are inhibited in 
a translationally invariant system [8]. This sequence of geometrical structure comes from the competition between 
the surface and Coulomb energies, indicating that these structural transitions are of first order. The mechanism of 
quark-hadron phase transitions at nonzero densities, however, has not yet been established by full QCD calculations, 
e.g., by Monte Carlo calculations of lattice gauge theory. 

Kinetics of the quark-hadron phase transitions has been studied primarily in the context of the early universe and 
ultrarelativistic collisions of heavy ions [9]. The astrophysical situation in which such kinetics could be significant 
is compression of dense stellar matter during collapse of a massive star’s core, spin-down of a neutron star, or mass 
accretion onto a neutron star from a companion star. The deconfinement transition in hadronic matter not containing 
strangeness, if proceeding via homogeneous nucleation [10], begins with appearance of a critical-size droplet of two- 
flavor (up and down) quark matter induced by fluctuations. Since a precritical droplet moves back and forth on a time 
scale for strong interactions, ~ 10“^^ s, which is many orders of magnitude smaller than that for weak interactions, 
flavor must be conserved during the precritical situation. As a consequence, a critical droplet of three-flavor quark 
matter containing strange quarks is unlikely to occur in spite of its stability over two-flavor quark matter. It is 
expected that a critical droplet of up and down quark matter appears either via thermal activation of the energy 
barrier separating the initial metastable (hadronic) phase from the stable (quark) phase in configuration space at high 
temperatures typical of matter created in highly energetic heavy-ion collisions and in stellar collapse, or via quantum 
penetration of this barrier at low temperatures appropriate to neutron star cores. We also note that during the growth 
of the critical droplet, weak processes producing strange quarks proceed fully in the stellar interiors in contrast to the 
case of heavy-ion collisions. 

In our recent papers [11,12], the time required to form a two-flavor quark matter droplet in a neutron star core, 
originally composed of nuclear matter in (3 equilibrium, was calculated for pressures near the point of the static 
deconfinement transition by using a quantum tunneling analysis incorporating the electrostatic energy and the effects 
of energy dissipation. It was found that the dissipation effects, governed by collisions of low-energy excitations in the 
hadronic phase with the droplet surface, significantly increase the degree of overpressure needed to form a droplet in 
the star, and that the formed droplet develops into bulk matter due to the screening of leptons on the droplet charge. 
We also discovered that the nucleation is likely to proceed via quantum tunneling rather than via thermal activation 
at temperatures below ~ 0.1 MeV typical of neutron star matter. It is of interest to note that strangeness may be 
present in the hadronic phase in the form of hyperons (A, etc.) [13,14] or in the form of a kaon condensate 

[15]. The presence of strange quarks in hyperons or kaons, altering not only the equation of state and the composition 
before and after deconfinement but also the quark-hadron interfacial energy, may affect the nature of the formation 
and growth of a quark matter droplet. 

A new physics content that we consider in this paper is the influence of the presence of hyperons on the dynamical 
deconfinement transition that may occur in the hadronic core of a neutron star accreting matter from a companion 
or rotating down. We particularly evaluate at what pressure a quark matter droplet forms via quantum tunneling in 
the metastable phase of hadronic matter containing hyperons (hereafter referred to as hyperonic matter). For this 
purpose, the thermodynamic properties of hyperonic matter such as the equation of state and the composition are 
estimated from a relativistic mean-field theory; we use the one with the parameters incorporating a binding energy of 
A in saturated nuclear matter [16]. By using an MIT bag model for the thermodynamic properties of quark matter, 
we calculate not only the pressure at which a deconfinement transition from hyperonic matter in /3 equilibrium takes 
place statically subject to flavor conservation, but also the corresponding jump of baryon density. In the overpressure 
regime, the time needed to form a quark matter droplet is obtained using a theory of quantum nucleation developed by 
Lifshitz and Kagan [17]; we build into this theory the effects of Coulomb energy, special relativity, energy dissipation 
[18], and dynamical [19] and static compressibility. Some of these effects are found to be crucial for the formation and 
growth of a quark matter droplet. We then estimate the critical overpressure required to form the first quark matter 
droplet in a neutron star core that consists of /3-stable hyperonic matter being compressed during a time scale for the 
spin-down or accretion. The resulting crossover temperature from the quantum tunneling to the thermal activation 
regime is shown to be large compared with the typical temperatures of matter in the core. By comparing the obtained 
critical pressure with the central pressure of the star with maximum mass, we evaluate the range of the bag-model 
parameters, i.e., the bag constant, the QCD fine structure constant, and the strange quark mass, where quark matter 
is expected to occur via deconfinement in the star. 

In Sec. II, the formalism to calculate the rate of nucleation of a stable phase in a metastable phase via quantum 
tunneling is discussed taking into account the effects of relativity, energy dissipation, and dynamical compressibility. 
In Sec. Ill, the static properties of the deconfinement transition are examined using a relativistic mean-field theory of 
hadronic matter with and without hyperons and a bag model for deconfined matter. In Sec. IV, on the basis of the 
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static properties obtained in Sec. Ill and the formalism described in Sec. II, the time needed to form a quark matter 
droplet in the metastable phase of nuclear or hyperonic matter is calculated in the overpressure regime; the critical 
overpressure and the quantum-thermal crossover temperature are also estimated allowing for the stellar conditions. 
Conclusions are given in Sec. V. In the Appendix, the quark distributions inside a droplet are considered. 


II. QUANTUM NUCLEATION THEORY 

In this section we first summarize a theory of quantum nucleation of a stable phase in a low-temperature first-order 
phase transition as advanced by Lifshitz and Kagan [17]. We then extend this theory to a relativistic regime where the 
effective mass of a critical-size droplet of the stable phase is comparable to the height of a potential barrier. Effects 
of energy dissipation [18] and dynamical compressibility [19] in the metastable phase, ignored above, are also built 
into the formalism to calculate the rate of quantum nucleation of the new phase. We assume in this section that the 
stable and metastable phases are electrically neutral quantum liquids with a single component. 


A. Lifshitz-Kagan theory 


Quantum tunneling nucleation of a stable phase in a first-order phase transition was first investigated by Lifshitz 
and Kagan [17]. Their analysis gives us a basic tool for evaluating the time needed to form a real droplet of the 
stable phase at low temperatures and at pressures in the vicinity of the phase equilibrium pressure. At such pressures, 
a nucleated droplet contains a large number of particles so that the resulting energy gain can compensate for the 
sharp energy increase in the interfacial layer. One may thus consider a droplet being in a virtual or real state to be 
a sphere macroscopically characterized by its radius R(t) and describe the tunneling behavior of a virtual droplet 
in the semiclassical approximation. By assuming that the velocity of sound Cs is sufficiently large compared with 
the velocity of the phase boundary, i.e., both phases are incompressible, the potential energy for a fluctuation of the 
radius R may be expressed in a standard form 

47rR3 

UiR)= g n2(A^2 - Ml) + . (I) 


Here is the chemical potential of the metastable (stable) phase calculated at fixed pressure P, n 2 is the 

number density of the stable phase, and Us is the surface tension. Expression (1) is derived from difference in the 
thermodynamic potential at fixed chemical potential between the initial metastable phase and the inhomogeneous 
phase containing a single droplet. As Fig. I illustrates, the potential barrier occurs between the initial metastable 
state and the state with a real droplet having a critical radius, Rc = 3as/n2{fJ.i — /r 2 ), which satisfies U{Rc) = 0. The 
height Uq of this barrier is given by Uq = ( 4 / 27 ) 47 rCTsi?g. 

The kinetic energy for a fluctuation of R is also necessary for the determination of the quantum nucleation rate. 
During such a fluctuation, the density discontinuity between the stable and metastable phases induces a hydrodynamic 
mass flow in the medium around the droplet. The velocity field v(r) is obtained from the continuity equation and the 
boundary condition at the droplet surface as 


v(r) 



r > R , 


0 , 


r < R , 


( 2 ) 


where R is the droplet growth rate, and ni is the number density of the metastable phase. The kinetic energy K{R) 
is thus given by 

KiR) = ^M{R)R^ , (3) 

where M{R) is the effective droplet mass. 


M (R) = Airpi 



( 4 ) 


Here pi is the mass density of the metastable phase. 
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At zero temperature, the time required to form a single droplet in the metastable phase may be calculated within 
the WKB approximation [20] by using the Lagrangian for the fluctuating droplet 

L{R, R) = ^M{R)if - U{R) . (5) 

Here no energy dissipation in the medium around the droplet is considered. The energy Eq for the zeroth bound state 
around i? = 0 is obtained from the Bohr quantization condition 

I{Eo) = , (6) 

where I{E) is the action for the zero-point oscillation 

fR- , _ 

I{E) =2 dR^/2M{R)[E -UiR)] , (7) 

Jo 

with the smaller classical turning radius R-. The corresponding oscillation frequency vq and probability of barrier 
penetration pq are given by 


dl 

dE 


E—Eq 


and 


Po = exp 


A{Eo) 

n 


( 8 ) 

(9) 


where A{E) is the action under the potential barrier 


A{E) = 2 



dRy^2MiR)[U{R) - E] , 


( 10 ) 


with the larger classical turning radius The formation time r is finally calculated as 

T = {i^oPo)~^ ■ 


( 11 ) 


B. Effect of relativity 

In case the potential barrier height Uq divided by is as large as the effective mass M{Rc) of a droplet of critical 
size, the effect of relativity on the quantum nucleation should be taken into account so that |i?| < c may be satisfied. 
Since this situation will appear in Sec. IV, it is instructive to describe a relativistic version of the Lifshitz-Kagan 
theory. Let us here assume that the velocity field itself is small enough to ensure the nonrelativistic description given 
by Eq. (2). The Lagrangian may then be rewritten as 


L{R, R) 


-M{R)& 


\ 


1 - 



-k M{Ry 


U{R) , 


( 12 ) 


where M{R) and U{R) are given by Eqs. (4) and (1), respectively. 

We proceed to obtain the time r needed to form a droplet from the Lagrangian (12) in the semiclassical approxi¬ 
mation. The Hamilton-Jacobi equation is derived from the Lagrangian (12) in a usual way as 


{Mc^f 





(13) 


where S{R,t) is the action associated with the Lagrangian (12). We perform the first quantization of Eq. (13) by 
replacing S with Ti/i^ and we thereby obtain a time-independent equation for the wave function ^’(1?) representing 
the state of energy E as 
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(14) 



(P 

+ {U- E){2Mc^ + E-U) 

dn^ 


^ = 0 . 


Since Eq. (14) bears a resemblance to the nonrelativistic Schrodinger equation, we can derive the semiclassical solution 
to the tunneling problem described in terms of Eq. (14) by following a usual line of argument developed in the 
nonrelativistic case [ 21 ]. As a consequence, a set of equations for determining the formation time r is obtained in the 
forms analogous to Eqs. ( 6 )-(ll). The energy Eq for the zeroth bound state around i? = 0 is now determined from 
the Bohr quantization condition 

I{Eq) =2tt (^Q + ^n . (15) 


Here I{E), the action for the zero-point oscillation, is rewritten as 

2 _ 

I{E) = - dR^y[2M{R)c‘^ + E- U{R)][E - U{R)] , 

c Jo 

with the smaller classical turning radius i?_; mo is the integer defined as 


Too 


IjEmin) 1 

2TTh 4 


(16) 


(17) 


where Ej^in is the maximum value of U{R) — 2M{R)c^, and [• • •] denotes the Gauss’ notation. Ey^in, being positive 
definite, yields the lower bound of the energy region where positive-energy states occur; by these we denote the states 
which ensure 2M{R)c^ -\- E — U{R) > 0 for arbitrary R. The zero-point oscillation frequency vq is determined by 
substituting Eq. (16) into Eq. ( 8 ). On the other hand, the probability po of barrier penetration ai E = Eq is given 
by Eq. (9) in which we now use the underbarrier action 

A(E) = - [ ^ dRj[2M{R)c'^ +E- U{R)][U{R) - E] , (18) 

c Jr_ 


with the larger classical turning radius R+. The obtained results for vq and po lead to the formation time r via Eq. 
(11). Equations (12)-(18) reduce to the nonrelativistic counterparts in the limit of c ^ oo. The higher order effect of 
relativity acts to raise the energy level and hence to enhance the tunneling probability, as will be clarified in Sec. IV. 


C. Effect of energy dissipation 


So far it was assumed that the system adjusts adiabatically to the fluctuation of R, i.e., a virtual droplet of the 
stable phase fluctuates in a reversible way. In a realistic situation, however, relaxation processes involved proceed 
at a finite rate; the density adjustment with varying R is necessarily accompanied by appearance of excitations in 
the metastable phase. These excitations give rise to dissipation of the total energy of the droplet until the system 
reaches a complete thermodynamic equilibrium. If the mean free path I of the excitations is much smaller than R, 
the excitations play a role in viscous transport of momentum from the high to the low velocity region with the rate 
of energy dissipation, 


dE 

dt 


Anri 


\ J’’' J r=R+0 


— Idirr] 



(19) 


where 77 is the viscosity of the metastable phase and v(r) is given by Eq. (2). This energy dissipation is thus built 
into the equation of motion for the droplet derived from Eq. (5) or (12) as an Ohmic friction force 

F = -167177 ^1 - RR . (20) 


In the case in which I ~ R, however, the hydrodynamic description of the dissipative processes ceases to be valid; 
considerations of the energy dissipation require the use of the corresponding quantum kinetic equation. In the regime 
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I ^ R where the excitations behave ballistically, they collide with the droplet surface and eventually relax in the 
medium far away from the droplet. The resulting momentum transfer dissipates the total energy of the droplet at a 
rate 


dt 


— IdTrarj 



( 21 ) 


where a is a factor of order unity that depends on the nature of the excitations and their interactions with the droplet 
surface, and thus exerts an Ohmic friction force on the droplet: 

F = -WTrar] (1 -. (22) 

V ni J I 

It is to be noted that expressions (21) and (22) are based on the hydrodynamic equation for the velocity field, which 
is not applicable to the situation considered here, and hence are not completely obvious. 

The influence of the Ohmic dissipation described above on the quantum nucleation rate was investigated by Bur¬ 
mistrov and Dubovskii [18] in terms of a path-integral formalism. This formalism was pioneered by Caldeira and 
Leggett [22] to consider the dissipation effects on quantum tunneling in macroscopic systems; it was noted that these 
effects act to reduce the tunneling probability exponentially. Hereafter, we extend the estimation of the nucleation 
rate to the system of finite temperature T with the help of path integrals, as advanced by Larkin and Ovchinnikov 
[23]. The probability p{T) of formation of a droplet in the medium undergoing the energy dissipation is then given 
by the expression similar to Eq. (9): 


p = exp 


MTV 

n 


(23) 


Here A{T) is the extremal value of the effective action ^eff specified in terms of an imaginary time t': 


S^siRir')] = J dr' |M(i?)c2^ 1 -b - M{R)c^ + U{R) 

: d ^ ff, ^ m 2 \ 

47r J-hp/2 sm 7r(T' - t")/TiI3 j 


(24) 


where /? = (ksT) ^ is the reciprocal temperature, Rt' denotes R{t'), and 7 (i?), the quantity determining the nonlocal 
dissipation term, is evaluated from the Ohmic friction forces (20) and (22) as 


7(^) = 


1 - 


r 8^1/2 

3 

2(7ra)^/^ 


m 




1-^ 

ni 


11/2 t I . 


(25) 


The effective action (24) is accompanied by the periodic boundary condition R{fif3/2) = R{—fif3/2). 

Generally, one obtains two types of trajectories that extremize the effective action SeS- One of these trajectories 
corresponds to a r'-independent classical trajectory R{t') = 2i?c/3, along which the action amounts to hUo/ksT. The 
nucleation described by this trajectory proceeds via thermal fluctuations. The other trajectory depends explicitly on 
t' and denotes the nucleation occurring via quantum-mechanical fluctuations. This statement is evident from the fact 
that the corresponding extremal action A{T) reduces to A(Eo) given by Eq. (18) in the limit of no energy dissipation 
and zero temperature. According as A{T) dominates over hUo/ksT or not, the nucleation of the stable phase proceeds 
via thermal activation or via quantum tunneling. There exists an abrupt transition from the quantum tunneling to 
the thermal activation regime at the temperature Tq derived from the relation 

.4m). ,26, 

Eor the purposes of practical application, it is useful to summarize the results for A(T) in the ballistic {R ^ 1) and 
nonrelativistic regime as given in Ref. [18]. In the case of weak dissipation, the dissipative processes add a relatively 
small quantity Ai{T) to A{Eo)-, the analytic expressions for A{Eo) and Ai{T) are given at Aq = 0 by 
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and 


A{Eo = 0) 


5-\/27r^ 

16 


1 - 


V2l 

ni 


y/asPiRl^^ 


Ai{T) = Alia 


n2 \ ^ yRj 

ni) I 


— + 0{T^) 

TT 


(27) 


(28) 


with J 2 « 1.49. For temperatures up to the crossover point Tq, not only does A{Eo = 0) retain its form (27) but also 
the term of order in Eq. (28) is negligible. For normal Fermi excitations in the medium, we obtain ij/l ~ nipp 
with the Fermi momentum pp = . Here it is instructive to note that in the case of weak dissipation, the 

Fermi velocity vp = nipp/pi is negligibly small compared with the characteristic velocity ^y2Uo/M{Rc) of a virtual 
droplet moving under the potential barrier. In the opposite case of strong dissipation [uf \/‘^Uo/M{Rc) for normal 
Fermi excitations], the dissipation effects control the dynamics of the droplet in configuration space and thus allow 
one to leave out the kinetic term in the effective action (24). The corresponding expression for A{T) reads 

AiT) = 4 ™ (^1 - ^s(r) , (29) 


where s(T), the normalized underbarrier action calculated along the extremal trajectory, is « 1.3 in the temperature 
range including T < Tq. 


D. Effect of dynamical compressibility 


In a system where the metastable phase is more or less compressible, a finite compressibility acts to reduce a 
portion of the liquid taking part in the hydrodynamic mass flow and hence the overall kinetic energy, leading to 
an exponential increase in the quantum nucleation rate. These effects were considered by Korshunov [19] in the 
nonrelativistic regime (c ji?j) by taking account of the time dependence of the velocity field. We formally extend 
his theory to the relativistic regime by replacing the kinetic term in the effective action (24) with 


Sm[R{t')] = 


rhpi2 


dr' 


i-hp/2 



X 



drATTpir"^ 




(30) 


where ip{r, t') is the velocity potential assumed to be spherically symmetric. The form of Sm adopted here is available 
when the flow velocity juj is much smaller than c^. The extremalization of Sm in P results in the wave equation 




0 

dr \ dr ) 


(31) 


Equation (31) can be solved under the boundary condition at the droplet surface based on the assumption that the 
stable phase is incompressible: 


ni 


- 1 i?= - 


dp{r, t') 


dr 


r—R{T') 


(32) 


It is obvious from Eq. (30) that Sm reduces to the kinetic term in Eq. (24) (hereafter defined as Sp) in the incom¬ 
pressible limit (ji?j <C Cs). The difference Sc = Sm — Sp has a nonlocal character similar to the dissipation term in 
Eq. (24), as shown in Ref. [19]. 

The compressibility term Sc, if the trajectory involved is taken to be the extremal one obtained at T = 0 and 
i?o = 0 in the nonrelativistic and dissipationless limit, leads to the expression [19] 


AciT) = -490 





Pi 


0's Cs 


keT 

n 


4 


(33) 
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Expression (33) describes just the principal term in the expansion with respect to T/Tq and yj2U[)IM{Rc)lcs. It 
is noteworthy that Ac{T) is proportional to {T/Tq)* for T ^ Tq. This behavior, attributable to the emission of 
the sound wave in the course of a fluctuation of the droplet radius, persists in the case in which the condition 
y/2Uo/M{R,) < c« is violated, the energy dissipation takes effect, and/or the droplet surface moves relativistically, 
as long as the time of flight of the extremal trajectory is finite for T — 0. Here we refrain from proving such finiteness 
explicitly, but this is a feature generally held by problems with a mass varying with R. 


III. STATIC DECONFINEMENT TRANSITION 

We now estimate at what pressure hadronic matter in (3 equilibrium undergoes a static deconfinement transition 
at zero temperature subject to flavor conservation.|^ For this purpose, we first describe the bulk properties of the 
/9-stable, zero-temperature hadronic matter with and without hyperons within the framework of a relativistic nuclear 
field theory [24]. Its Lagrangian density, denoting the interactions between baryons by the exchange of A4 mesons 
(A4 = a,uj,p), is given by 

-rUB + gaBO- - - l^gpBltiT • P^ji’B 

-U{a)- mi)'ipi , (34) 

I 


with 


tOpy = - d^ujfj, , 

PfJ,U ~ dpPv ~ dj^Pp ■ 

Here '<pB{i) is the Dirac spinor for baryons B = p, n. A, E+, S“, E°, S“, (for leptons I = t is the isospin 

operator, gMB is the coupling constant between B baryons and Ai mesons, and rrii is the rest mass of the particle 
of species i {i = B,A4,l). The potential U{a), denoting the self-interactions of the scalar field and playing a role in 
reproducing the empirical nuclear incompressibility, takes a specific form, 

U{a) = ^airrinigcrNCff + ^ 02 ( 50 -at o')'^ , (35) 

where N represents the nucleons {N = n,p). 

From the Lagrangian density (34) we obtain a couple of the Euler-Lagrange equations with respect to the meson 
and baryon fields. In the mean field approximation adopted here, the meson fields are replaced by their expectation 
values and are determined from the classical field equations [14] 

mlujo = ^ guiBUB , (36) 

B 


^Ip3o = 

B 




(37) 


mla = - 


dU (tr) 
da 


1 

27r^ 


B 


tB\J 


1 -b t R — In 




(38) 


where is the number density of B baryons, Mb = mB — g^Bf^ is the effective mass of the baryon, t^b is the 
third component of isospin of the baryon, and ts = kp^B/MB with the corresponding Fermi wave number kp^B = 


^Hereafter, we take units in which % = c = kp ~ 1- 
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(STT^ns)^^^- Here the meson fields are constant for any space and time since the considered system is static and 
uniform. This system is also isotropic, leading to u; = Pa = 0. We have assumed, furthermore, that the charged 
components of p mesons as well as the pions and the kaons have vanishing expectation values, as is normally the case. 

The total energy density £tot and pressure P for the hadronic system are then calculated from the energy-momentum 
tensor derived from the Lagrangian density (34) as a sum of the hadronic term and the leptonic term [14]. The resulting 
expression for £tot is 


with 


£tot — £h + £l , 


Eh (a) + 




{2t% + 


1 -I- t R — In 


ts 



(39) 


(40) 




(2ti + l)t/-y/l t‘f — In \ ti y 1 1 


(41) 


where ti = kF,i/Tni with the Fermi wave number of I leptons kp^i = (37r^ni)^/^. The pressure is obtained as 


P = Ph + Pl , 


(42) 


with 


Ph 


B 




(43) 


Pl = -El + ^ nipi , 

i 


(44) 


where 





+ + ffwBWo -I- gpBLsBPSO 


is the chemical potential of B baryons, and 


(45) 


gi = 



(46) 


is the chemical potential of I leptons. Here electrons and muons are looked upon as ideal gases; the density region 
considered here is high enough for the kinetic energy to dominate over the energy induced by Coulomb interactions. 
The pressure P is in turn related to the total energy density £tot, equivalent to the mass density p, via 




(47) 


where nt = ub is the total baryon density of hadronic matter. 

The parameters g^B, gpB, guiB, oi, and 02 contained in Eqs. (34) and (35) have been taken from the values 
determined by Glendenning and Moszkowski [16]. The properties of saturated nuclear matter adopted in Ref. [16] 
are characterized by the binding energy per nucleon (16.3 MeV), the saturation density (0.153 fm“^), the symmetry 
energy (32.5 MeV), the incompressibility (300 MeV), and the effective nucleon mass (0.7m„). In such a way as to 
reproduce these properties, the nucleon-meson coupling constants and the parameters determining the strength of the 
a self-interactions were chosen as PcrAr/wcr = 3.434 fm, g^jN/m^ = 2.674 fm, PpAr/wp = 2.100 fm, oi = 0.00295, and 
02 = —0.00107. The hyperon-meson couplings, denoted by the ratios 


— 


gaH 

g^N 




guiH 

guN 


^pH 


gpH 

gpN 


(48) 
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were determined mainly from the empirical A binding energy (28 MeV) in saturated nuclear matter as x^h = XpH = 0.6 
and = 0.653. Here the coupling constant ratios are assumed to be the same for all hyperon species. It was noted 
that the ratios thus chosen are compatible with the observational lower bound of the maximum neutron star mass 
and with the upper bound of x^h stemming from the fit to hypernuclear levels. 

In the determination of the equilibrium composition of the hadronic system at a given Uh, we need not only Eqs. 
(36)-(38), but also additional constraints. These are the charge neutrality condition 

Qbtib - ne - = 0 , (49) 

B 


where qb is the electric charge of the baryon of species B, and the /3-equilibrium conditions 


Mb = Mn - gsMe , (50) 

Mm = Me • (51) 


We have thus calculated the fraction Yi = nijnb of i particles {i = B,l) in matter with and without hyperons; the 
results have been plotted as a function of nt in Fig. 2. We can observe in this figure that the lepton fractions in 
hyperonic matter are drastically reduced by the appearance of hyperons in the density region nb ^ 2no. In Fig. 
3 we have depicted the equations of state of hadronic matter as calculated allowing for and ignoring the presence of 
hyperons. For comparison we have also included the equations of state of equilibrium nuclear matter in the many- 
body calculations [25,26] based on realistic nucleon-nucleon interactions. We thus see that the presence of hyperons 
at high densities as predicted from the relativistic mean-field theory softens the equation of state considerably. The 
relativistic field theory used here automatically ensures the causality condition that the sound velocity 


/^\ 1/2 ^ / nt, dP \ 
\dp J \p + Pdnb) 


(52) 


be confined within the velocity of light. 

We turn to the description of the bulk properties of uniform quark matter, deconfined from the /3-stable hadronic 
matter mentioned above, by using an MIT bag model. We begin with the thermodynamic potential of q quarks, where 
q = u,d, and s denote up, down, and strange quarks, expressed as a sum of the kinetic term and the one-gluon-exchange 
term [2,27,28] 


4^2 


dqidq - + ^m^ln 


Mg 


+ (m, - to2)1/2 



- ml In 

- 2(^2 -^ 2 ) 2 } ^ 




(53) 


where and Mg are the q quark rest mass and chemical potential, respectively, and as denotes the QCD fine structure 
constant. For u and d quarks, the ratios of the mass to the chemical potential are negligibly small; we thus use the 
following expressions in the massless limit: 


flit 




fid 


\ TT 


(54) 


The number density Uq of q quarks is in turn related to fig via 


riq = 


9flg 

dpq 


Then, the total energy density for the quark system is 


£tot — ) 


(55) 


(56) 


with 
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£Q — 'y 'X ^q l^q'^q) 

9 

2a, 


= 1 - 


TT / 47r2 


‘5^4 

+ 1) - ln(x, + r?,)] 


asirq 

2tt^ 


{2xl{xl + 2?7^) - 3[xsr]s + ln(xs + rj^)] } + b , 


(57) 


where Xg = — rngjiris, rjg = %/!■ + a;2, is given by Eq. (41) and b is the energy density difference between the 

perturbative vacuum and the true vacuum, i.e., the bag constant. The total baryon density is now given by 


q 


By using the thermodynamic relation, we obtain the pressure 


P = Pq + Pl, 


with 


Pq = - - b , 

q 


(58) 


(59) 

(60) 


where Pl is given by Eq. (44). 

The parameters mg, ag, and b are basically obtained from the fits to light-hadron spectra. However, the values 
of these parameters yielded by such fits do not correspond directly to the values appropriate to bulk quark matter 
[28]. Taking account of possible uncertainties, we set 50 MeV fm“^ < b < 200 MeV fm“^, 0 < nig < 300 MeV, and 
0 < tts < 1. The ranges of ag and mg so chosen are consistent with the results of Barnett et al. [29] that are inferred 
from experiments and renormalized at an energy scale ~ 1 GeV of interest here. 

In obtaining the composition of the quark system, we take note of its relation via flavor conservation with that of 
the /3-stable hadronic system [30]: 


E„ \ 

1 / 2 

1 
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Yd 
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\y:-) 


(61) 


where Yq = Uq/rii, with nj, given by Eq. (58) is the fraction of q quarks in the deconfined phase. The absence of leptonic 
weak processes and the charge neutrality condition ensure that the fraction Yi of I leptons remains unchanged before 
and after deconfinement. 

Let us now proceed to obtain the pressure Pq of the static deconfinement transition from the relation 


/4fc (Tq) [hadronic phase — (Eg) [deconfined pha 


(62) 


where 


Mb 


£tot + P 
Hb 


(63) 


with the corresponding total energy density £tot and total baryon density Ub is the baryon chemical potential at fixed 
P. We have calculated the pressures Pq for the parametric combinations between 50 MeV fm“^ <b < 200 MeV fm“^, 
rUg = 0,150,300 MeV, and 0 < < 1, allowing for and ignoring the presence of hyperons in the hadronic phase. 

Eigure 4 illustrates the resultant contour plots of the pressures Pg on the b versus ag plane. As specific values of Pg 
shown in Fig. 4, we have chosen the pressure P^- at which hyperons appear in the hadronic phase, the central 
pressure Pmax of the maximum-mass (Mmax) neutron star having a hadronic matter core, and the central pressure 
Pi ,4 of the star (not including quark matter) with canonical mass 1.4 Mq; these values are tabulated in Table I. P^- 
can be determined from the relation coming from Eq. (50): 
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=0 Me 5 


(64) 


Mn — Me- \n^- 

where In^,. =o denotes the energy of the lowest state for a 11“ hyperon. The tabulated values of Pmax and Pi .4 have 
been calculated from the structure of the star whose core consists of /3-stable nuclear matter or hyperonic matter. In 
these calculations, the gravitational mass of the star observed by a distant spectator has been obtained as a function 
of the central pressure or mass density of the star from the general-relativistic equation of hydrostatic balance in the 
nonrotating conhguration, i.e., the Tolman-Oppenheimer-Volkoff equation [31], and from the corresponding equation 
of state. The equation of state used here has been extrapolated to lower densities appropriate to matter in the 
crust made up of a lattice of nuclei embedded in a sea of electrons and, when present, in a sea of neutrons. Such 
extrapolations make only a negligible difference in determining the stellar structure, since the crust contains mass of 
^ O.OIM 0 , confined within a few percent of the total mass [32]. 

The effect of a finite s quark mass acts to destabilize the quark system and hence to lower the contours of various 
Pq on the b versus as plane, as shown in Figs. 4(b)-4(d). We have found that for rus ^ 200 MeV, such lowering is 
fairly small compared with the adopted range of b and as- Thus, we shall generally take TOs to be zero in describing 
the bulk properties of quark matter. We also observe in Figs. 4(a) and 4(b) that the contours of Pq = Pmax and 
Pq = Pi ,4 reveal only a weak dependence on the presence of hyperons. This result implies that within the confines of 
the present models for matter, the degree of strangeness-induced stabilization of the hadronic phase compensates for 
that of the deconfined phase as long as rUs ^ 200 MeV. 

For the parameters as = 0.4, b = 100 MeV fm“^, and iris = 0, we have depicted in Fig. 5 the baryon chemical 
potentials and total energy densities for three electrically neutral bulk phases. The first phase consists of /3-stable 
hadronic matter with or without hyperons. The second phase is composed of the matter deconfined therefrom 
according to Eq. (61). Lastly, for comparison, we consider a phase of u,d, and s quark matter in /3 equilibrium: 
IJ-d = fJ-uP fJ'e and fj,^ = fds- Hereafter, we shall refer to this /3-stable quark matter as strange matter.^ It is instructive 
to note that strange matter composed of massless quarks satisfies = 12,1 = ris and rie = = 0. The upper crossing 

shown in Fig. 5(a) [5(b)] denotes the deconfinement pressure Pq, at which the total baryon density changes as can 
be seen from the Maxwell equal-area construction plotted in Fig. 5(c) [5(d)]. The density discontinuity for hyperonic 
matter is significantly smaller than that for nuclear matter, a feature with immediate relevance to the nucleation 
of the deconfined phase as will be discussed in Sec. IV. We likewise observe that the presence of hyperons in the 
hadronic phase suppresses the chemical potential difference between deconfined and strange matter. This is because 
the dominance of A hyperons in hyperonic matter renders the quark fraction Yq in deconfined matter similar to that 
in strange matter. 

We conclude this section by indicating that thermal effects on the static properties of the deconfinement transition 
may be safely omitted at low temperatures appropriate to neutron star cores. The zero-temperature models for 
hadronic and quark matter, described above, hold over an even wider range of temperature, T < 10 MeV. This is 
because the relative chemical potentials — ^i\ni=o of i particles (/ = except for minor components with 

particle fraction Yi < 0.1 {i = P) or < 0.001 {i = l,q), take on a value of ~ 30-600 MeV, sufficiently large compared 
with the temperature. This temperature range will thus be considered below. 


IV. DYNAMICAL DECONFINEMENT TRANSITION 

In this section we first construct the effective action for the fluctuational development of a quark matter droplet in 
the metastable phase of hadronic matter with and without hyperons. We thereby estimate the time necessary for the 
quantum-mechanical formation of the first droplet in a neutron star as a function of pressure. The critical overpressure 
is determined in such a way as to make this formation time comparable to the time scale for the compression due to 
the star’s spin-down or accretion; at this overpressure, the quantum-thermal crossover temperature is estimated. We 
finally determine the window of the bag-model parameters &, as, and rUs where quark matter is expected to nucleate 
in the star. 


^ Although strange matter is not related to the deconfinement transition of interest here, it would occur once a droplet of the 
deconfined matter forms and grows into bulk matter in a neutron star. This is because the weak process u + d —> u + s converts 
the bulk deconfined matter into strange matter during the growth of the droplet [1]. 
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A. Effective action 


The quantum nucleation theory described in Sec. II is directly applicable to the case in which the metastable and 
stable phases are electrically neutral quantum liquids with a single component. The deconfinement phase transition 
of interest here, however, involves matter with multiple components such as baryons, leptons, and quarks. Moreover, 
the inhomogeneous state of matter including a single quark matter droplet leads inevitably to a violation of the local 
charge neutrality, a property that is inherent in the systems (e.g., nuclei) where strong interactions compete with 
Coulomb interactions. It is thus necessary to generalize the quantum nucleation theory by taking into account the 
multiplicity of components and the electrostatic energy. 


1. Nuclear matter 


Let us now construct the effective action, leading to the rate of quantum nucleation of quark matter in the case 
in which the metastable phase consists of /3-stable nuclear matter, by following a line of argument of Refs. [11,12]. 
We begin with the potential energy U{R) which is defined in the limit of Cg 3> |i?| as the minimum work needed to 
form a quark matter droplet of radius R in the metastable phase. Pressure equilibrium between quarks and nucleons 
embedded in a roughly uniform sea of leptons induces the excess positive charge inside the droplet, because the baryon 
density of the quark phase is larger than that of the hadronic phase by ~0.1-0.7 fm“^. The charged components are in 
turn distributed in such a way as to screen the droplet charge. The screening efficiency of i particles (i = e, ii,p, u, d) 
may be partially measured from their Thomas-Fermi screening length 


Atf.i 


1 

dru ) 


(65) 


where Qi is the electric charge of the particle of species i. For the droplet sizes of interest, R ^ 10 fm, we can 
assume the fluid of protons to be homogeneous since Atf,p ^ 10 fm in the considered density range nt ^ 4no. The 
other components, whose screening lengths are confined within 10 fm, deviate more or less from uniformity. On the 
other hand, the flavor conservation, holding Hxed the ratio of the number of u quarks to that of d quarks inside the 
droplet, makes the screening action by quarks ineffective, since u and d quarks have electric charge with opposite 
sign. We assume that u and d quarks are distributed uniformly; it turns out that the quark screening, i.e., the static 
compressibility of the fluid of quarks, does not contribute signihcantly to a reduction of the potential energy of the 
droplet (see the Appendix). 

For the initial metastable phase of nuclear matter in /3 equilibrium under pressure P, we can obtain the total energy 
density etot,H, the baryon density rib^H, and the baryon chemical potential pb.H by using Eqs. (39), (42), and (63). 
We then express the energy density for the inhomogeneous phase containing a single droplet as the sum of bulk, 
interfacial, and Coulomb terms: 


£d(f) = 9{R - r)eQ{nu,nd,ns = 0) -I- 9{r - i?)£pp(n„, Up, ub^n = 0) 

+£L[ni{r)] + es{r) + ^E(r)^ . (66) 

OTT 

Here r is the distance from the center of the droplet, eg, Eh, and el are given by Eqs. (57), (40), and (41), respectively, 
es{t) is the increase in energy density due to the quark and nucleon distributions in the interfacial layer whose thickness 
we assume to be much smaller than i?, and E('r) is the electric field. The baryon density for the inhomogeneous phase 
is given by 


nb,D{r) = 9{R - r)nb,Q + 9{r - R)nb,H , (67) 

where rib^q is the baryon density inside the droplet determined by Eq. (58). We thus obtain the potential energy 
by integrating the difference in the thermodynamic potential per unit volume at chemical potential ttb,H between 
the initial metastable phase and the inhomogeneous phase over the system volume V, which is related to P via 
P = -d{VEtot,H)/dV, as 


U{R) 


/ dV{EDir) - Etot,H - £'(’’) - nb,H]} 

Jv 


IV 
AttR^ 


nb,Qitlb,Q — Pb,H) + /S.El{R) + ‘I'KagR^ + Ec{R) 


( 68 ) 
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where 


M 6 ,Q 


£Q{nu, rid, rzg = 0) + eL{ni,H) + P 
nb,Q 


is the chemical potential for the quark phase and 


/\El{R) = [ dV{eL[ni{r)] 
Jv 


(69) 


(70) 


is the excess of the lepton energy over the initial value. Here tii^h is the number density of I leptons in the initial 
metastable phase of nuclear matter, is the surface tension, and Eq is the electrostatic energy. 

Since the quark-hadron interfacial properties are poorly known at finite densities, the expression for ag has been 
taken from the Fermi-gas model for a quark matter droplet in vacuum [33] as 


CTs 


3 J ’^q) 

47r I 6 

9 k 



1 

37r 




-2^qmq{nl ~ + ”^9 


^9 + (^9 “ 


nia 


(71) 


where is set to be zero. This expression arises from the reduction in quark density of states due to the presence of the 
droplet surface. The typical value of ag for matter without strangeness is obtained as cTg « (3/47r^)(TO„^^-|-md/r^) ~ 10 
MeV fm“^, since rUu ~ md ~ 10 MeV and ~ ~ 500 MeV. The foregoing Fermi-gas description includes neither 

0{ag) corrections nor curvature corrections. The former have yet to be determined, while the latter destabilize a 
quark matter droplet in vacuum by producing an energy 

i?curv = -^/r^ (72) 

TT ^' 

9 


for fj,q ^ niq and = 0 [34]. Instead of explicitly including these curvature corrections in the potential (68), we set 
the range of ag as 10 MeV fm“^ < CTs < 50 MeV fm“^ by incorporating into ag the increase in the interfacial energy 
yielded by the curvature energy (72) of the droplet having a radius in excess of ~ 3 fm. Shell effects as encountered in 
nuclei are likely to appear remarkably in light quark matter droplets having a baryon number < 10 [28], which cannot 
be well described by the Fermi-gas model for finite quark matter as adopted here. Nevertheless, we may assume this 
model to be useful, since the baryon number contained in a virtual droplet moving under the potential barrier is ^ 10 
(typically of order 100). The contribution of nucleons to the interfacial energy, ignored here, is expected to be much 
smaller than that of quarks at rib^H ~ no, but it is uncertain in the considered density range Ub^H ~ 4-6no. 

In the absence of the lepton screening on the droplet charge, 


Zoe 


AttR^Pqo 

3 


with 


PQo = e 


2tIu nd ng 

3 



we obtain Ec = and 


AEl = Zop,e,H , 


(73) 


(74) 


(75) 


where Pe,H is the chemical potential of electrons in the initial metastable phase of nuclear matter; typically, Pe,H ~ 
250-300 MeV. Here Eq. (75) reflects the fact that the global charge neutrality in the system volume V is guaranteed by 
the lepton gas which satisfies the /^-equilibrium condition (51). The excess lepton energy AEj^, which is proportional 
to R^, is naturally absorbed into the chemical-potential difference term in the potential (68). The number densities 
n„ and Ud are then determined from the flavor conservation (61) and the pressure equilibrium between quark and 
nuclear matter 


Pq=Ph 


AEl 
AttR^ ’ 


(76) 
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where Pq and Ph are given by Eqs. (60) and (43), respectively. In Eq. (76) the surface and Coulomb pressures arising 
from Ai^asR^ and Ec{R) have been ignored since these pressures are much smaller than the bulk pressures Pq and 
Ph-^ n-u and rid are thus independent of R. It is instructive to note that the flavor conservation (61) and the charge 
neutrality (49) in the hadronic phase result in the proportionality of pqo to rie^H + and to — 1: 

PQO = e{ne,H + ^ - 1) , (77) 

\nb,H J 


where the ratio rid^Q/rih^H amounts to 1.2-1.8. 

We next consider the screening action by leptons on the excess droplet charge Z^e. The number density n;(r) of 
I leptons, whose deviation Sni{r) from npH is small enough to be treated as a perturbation as will be shown below, 
has been determined from the linear Thomas-Fermi screening theory [35]. This determination [11] leads to 


Sni{r) 


dni 

dp-i 


z(j){r) = 




47re^ 


z(j){r) 


(78) 


where 4){r) is the electrostatic potential derived from the Poisson equation V‘^(j)(r) — = —47rpQo^(.R ~ ?') with 

the reciprocal of the screening length kl = J^ Us solution reads 


ArrpQQ 


(j){r) = < 


1 — exp(—Kii?)(l + klR) 


sinh(Kir) 


klt 


r < R , 


^^^[KLi?cosh(Kii?) - sinh(KLi?)]^^^^^-!^ , r>R. 

Kj HiLT 


The resulting electrostatic energy is given by 


Ec{R) = 


Pqo 


{-3 + [klRY + exp(-2KLi?)[3 + QklR + ^{ulRY + 2(kl7?)^]} . 


(79) 


(80) 


We have confirmed that the linear Thomas-Fermi approximation may be safely used; the validating condition Pe,H ^ 
e(j){r) is well satisfied, since e4>{r)/p,e,H ^ 0.1 for i? < 10 fm.^ 

Figure 6(a) exhibits the potential for the formation and growth of a quark matter droplet, characterized by 6 = 100 
MeV fm“^, Us = 0.4, and as = 30 MeV fm“^, in the initial metastable phase of nuclear matter at pressures above 
Po ~ 453 MeV fm“^. This figure contains three cases in which the electrostatic energy is set to be zero, the nonscreened 
value 3^0 6^/57?, and the lepton-screened value given by Eq. (80). By comparing these cases, we observe that the 
lepton-screened case shows no dominance of the Coulomb energy over U(R) for large R in contrast to the nonscreened 
case, and agrees fairly well with the case in which Ec = 0. This is because the cloud of leptons spreads inside the 
droplet and from its surface outward over a scale of the screening length ~ 5 fm (see Fig. 9 in the Appendix). 
We have thus confirmed a salient feature found in Ref. [11] that the lepton screening effects prevent the formation of 
a Coulomb barrier which causes the droplet to remain finite. The good agreement between the lepton-screened and 
Ec — 0 cases allows us to calculate the nucleation time without including Ec in the potential (68). Hereafter, we 
shall leave out Ec in the potential (68); the resulting potential is characterized by the critical droplet radius 


Rc 


_3o^_ 

nb,Q{p,b,H - Pb,Q) - PQ0Pe,H/e 


(81) 


and by the barrier height Uq = (4/27)47rtTsi?^. 

The effective mass M{R) of a quark matter droplet in /3-stable nuclear matter may be estimated from the kinetic 
energies of nucleons and leptons, in a way analogous to that described in Sec. II A. We obtain the nucleon kinetic 
energy of the form (3) by deriving the baryon velocity field from the baryon continuity equation and the 

boundary condition at the droplet surface; the expression for Kh is 


^As long as ct® oc the surface pressure vanishes. 

'^The energy increment yielded by the gradient of 5ni, present over a scale of the screening length ~ 5 fm inward and 
outward from the droplet surface, is neglected in the present study. 
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Kh = 27re//(n„, Up, ub^n) 
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(82) 



nb,Q 

nb,H 




By substituting the lepton distribution (78) into the lepton continuity equation, the lepton velocity field can be 
determined as 


ni^ne \kl J 


R 


ni^ne \kl J {kltY 

Ki,i?exp(—Kii?)[KLr'cosh(KLr) — sinh(KLr)] , 
(klR)'^ — KLi?sinh(Ki_R)(l + K^r) exp(—K^r) 
The lepton kinetic energy Kl, which is evaluated from Eq. (83) as 


Kl < 2TT£L{niM) 


ni 

. I ^ 


PQO 

ni,He 


R^R^ 


r < R , 
r > R . 


(83) 


(84) 


proves negligible, partly because eff(n„, rip, = 0) ~ 1000 MeV fm“^ is much larger than £L(ni,H) of order 

10 MeV fm“^ and partly because Eq. (77) and ne,H ^ '^h,h ensure {nb^Q/ub^H — 1)^ (^f/^i)(PQo/^i.ffe)^. The 

effective mass M{R) is finally obtained as 


M{R) = 4:TTeH{nn,np,nB^N = 0) 


nb,H J 


(85) 


The effect of energy dissipation discussed in Sec. II C plays a crucial role in determining the time needed to form 
a quark matter droplet in /3-equilibrated nuclear matter, as shown in Ref. [12]. Let us now assume that the nucleons 
and the leptons are in a normal fluid state. Then, elementary excitations of particle species i {i = N,l) appear as 
quasiparticles in the vicinity of their respective Eermi surfaces during a fluctuation of R. Since the mean free path 
of i excitations is estimated^ as k ^ — ^i|ni=o)^T“^ > 10^ fm, we can observe that each k is sufficiently 

large compared with the critical droplet radius Rc ranging typically 1-5 fm. Consequently, the nucleon and lepton 
excitations behave ballistically and collide with the droplet surface. The resulting transfer of momentum flux exerts 
on the droplet an Ohmic friction force of the form analogous to Eq. (22); this is written as a sum of nucleonic and 
leptonic terms 


with 


F = FhIub^n^O + Fl , 



( 86 ) 


(87) 


Fl = -WnfiKBR) 
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PQO 
Kl \ni,He 


aini^HkF,i,H 


R^R 


( 88 ) 


Here kF,i,H — is the Fermi momentum of I leptons in the metastable phase, ai is a factor of order unity 

that depends on the properties of i excitations and on their interactions with the droplet surface, and /{klR) is a 


® This estimate of U is based on the Landau theory of the Fermi liquids [36]. Here we have assumed that the coefficient 
Ci affixed to — M*ki=o)^T“^ is of order unity. However, Ci depends on the interaction between quasiparticles (i.e., 

the Landau parameters). For nucleons, such an interaction is repulsive in the high density region of interest here [37,38]; the 
values of Cjv range ~ 0.1-1. For leptons, being relativistic and degenerate, polarization effects of the medium on the exchanged 
photons [27,39] induce attraction between quasiparticles; the values of Ci are roughly 10. These considerations of Ci do not 
change the conclusion U Rc. 
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function that stems from the velocity field (83) occurring due to the lepton screening effects, behaves as oc (klR)^ 
for klR <C 1, and monotonically approaches unity as klR increases.^ Since ~ 250-450 MeV fm“^ is 

one order of magnitude larger than ~ 30-70 MeV fm“^, we ignore the leptonic contribution (88) to the 

Ohmic friction force. The Fermi velocity VF,H\nB^N=o averaged over nucleon species, where vf,h is defined as 


VF,H 


£H{nn,np,nB^N) 


(89) 


is ~ 0.4c, and hence VF.H\nB^N=o is in excess of y^2Uo/M{Rc) ~ 0.1-0.3c. As a consequence, by recalling how to 
distinguish between weak and strong dissipation as described in Sec. II C, we find that except for a^r <C 1, the 
nucleation time is determined by the dissipative processes rather than by the reversible droplet motion. 

The effects of relativity and dynamical compressibility as considered in Secs. II B and II D make only a negligible 
change in the estimates of the time necessary for the quantum nucleation of quark matter in the metastable phase 
of nuclear matter. This is because the droplet surface moves under the potential barrier at a velocity [typically 
~ ^y2Uo/M{Rc)] fairly low compared with c^ ~ 0.8c (nuclear matter) and Cg ~ 0.6c (deconfined matter). These 
comparisons of velocity scales lead us to the effective action of a form similar to Eq. (24): 



U{R) + a 1 



X 







sin^ 7r(T' — t")I[3 


(90) 


Here U{R) is taken from Eq. (68) in which we set Ec = 0, and an = o:p = a is assumed for simplicity. The extremal 
action A{T) for the quantum nucleation is obtained in a form identical to Eq. (29) as 

nb,Hj y 

where Rc is given by Eq. (81). 


A(T) = ‘iira 
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nNkF,N Rts{T) 


(91) 


2. Hyperonic matter 

Let us now construct the effective action for the development of a quark matter droplet in the initial metastable 
phase of /3-equilibrated hadronic matter with hyperons. In this construction, we may utilize expressions (65)-(89) 
developed for matter without hyperons by identifying the quantities characterizing the metastable phase at a given 
P, i.e., £tot,H, and ui^h, with those for matter with hyperons as well as by removing the constraints 

ub^n = 0 and Ug = 0. Notice that the lepton fractions plotted in Fig. 2(b) are sufficiently small at baryon densities 
of interest, rib^H 4-6no, to confine the ratio \pQo/nb^Qe\ of the electric charge to the baryon number inside the 
droplet within 10“^. This behavior stems from the fact that the ratio |pQo/''^b,Qe| is proportional to the sum of the 
lepton fractions {ne,H + np^H)/nb,H, as can be seen from Eq. (77). When using expressions (65)-(89), therefore, 
we omit the quantities associated with electricity by setting pgo = 0- Correspondingly, the screening action by the 
charged components is neglected. Such an omission is valid for the droplet sizes of interest here, R ^ 20 fm. We 
have confirmed that for R ^ 20 fm, the Coulomb energy continues to be trivial due mainly to the screening action 
by leptons and due partly to that by charged baryons. 

By noting Eqs. (68) and (69), the potential U{R) for the formation and growth of a quark matter droplet in /3-stable 
hyperonic matter has been written as 

47ri?3 

U{R) = —^nb,Q{pb,Q - Pb,H) + ^nagR^ , (92) 


®Note that Eq. (86) reduces to Eq. (8) in Ref. [12] in the absence of muons and in the limit of klR —>■ oo. Thus, the 
contribution of the electron quasiparticles to the friction force was rather overestimated in Ref. [12]. 
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with 


eqinq) + eL{ni,H) + P 

/ib,Q = - , (yoj 

nb,Q 

where nt^q, eg, and £l are given by Eqs. (58), (57), and (41). Here as, the surface tension, has been assumed to range 
from 10 to 50 MeV fm“^ by using the Fermi-gas expression (71); this range stems primarily from the uncertainty in 
TOs- The curvature term (72) contributes only a little to the interfacial energy of the droplet with typical R ^ 5-15 
fm. The quark number densities nu, Ud, and rig have been obtained from the flavor conservation (61) and the pressure 
equilibrium between quarks and baryons Pq = Ph, where Pq and Ph are given by Eqs. (60) and (43), respectively; 
modification of such pressure equilibrium by the interfacial energy is negligibly small for typical R. In Fig. 6(b) we 
have depicted the potential energy (92), calculated for b = 100 MeV fm“^, ag = 0.4, rUg = 0, Cg = 30 MeV fm“^, 
and various pressures above Pq « 163 MeV fm“^. The form of this potential is identical to the standard one denoted 
by Eq. (1); the critical droplet radius is given by 


Rc 


_Sog_ 

«b,Q(Mb,i? - Mft.g) 


(94) 


from which the potential barrier height Uq is determined as Uq = (4/27)47r(Tsi?^. 

The effective mass M(R) of a quark matter droplet in the metastable phase of hyperonic matter is then obtained 
from the baryon kinetic energy Kh given by Eq. (82) as 


M{R) = 47r£fl'(nB) 


nb,H ) 


(95) 


In this determination we ignore the lepton and quark kinetic energies since we have set pqo = 0. Expression (95) does 
not include the relativistic effect on the hydrodynamic mass flow of baryons. As will be shown in the next paragraph, 
the baryon velocity field, given by Eq. (2) in which ni = nb^H and n 2 = Ub^q, is far smaller than c, although the 
droplet motion itself is relativistic in many cases. 

In constructing the effective action for the system having a quark matter droplet in /3-equilibrated hyperonic 
matter, it is instructive to compare the velocities |i?|, Cg, and vf,h characterizing the droplet motion under the 
potential barrier. It is of great importance to note that over a wide range of the parameters Ug, b, as, and rUg 
with relatively large as (^ 0.2) and as {P 30 MeV fm“^), the density discontinuity Ub^q — Ub^n is small enough to 
make ^J2Uq/M{Rc) higher than 0.3c for typical Rc ^ 5-15 fm. We should, therefore, take into account the effect 
of relativity on the quantum nucleation as mentioned in Sec. II B. The velocity of sound, Cg ~ 0.5c, of hyperonic 
matter cannot fully surpass the rate \R\ of growth of a virtual droplet ranging typically ~ 0.1-lc. Thus, the effect of 
the dynamical compressibility of the metastable phase as discussed in Sec. II D is likely to have consequence to the 
quantum nucleation at temperatures near the quantum-thermal crossover point .[^ The absolute value of the baryon 
velocity field, which is less than |(1 — nb^q/nb^H)R\ [see Eq. (2)], is limited within 0.1c. This leads us to use the 
effective mass M{R) given by Eq. (95). The effect of energy dissipation is controlled by low-lying excitations of 
baryons induced by the fluctuation of R. These excitations, whose mean free paths are sufficiently large compared 
with typical Rc (except for minor components), collide with the droplet surface and exert the Ohmic friction force 
Fh given by Eq. (87) on the droplet. Here we have assumed the baryons to be in a normal fluid state, and we have 
ignored the negligible contributions of leptons and quarks to the energy dissipation. The Fermi velocity given 
by Eq. (89), averaged over baryon species, is estimated to be ~ 0.3c, so that \R\ > vf.h is satisfied in the relativistic 
regime. We thus confine ourselves to the dissipationless cases. We finally write the effective action as 


Ses[R{T')] = dr' \ (-1 + Vl + rA J 

J-/3/2 I ^ 'Ft 


1 


R(t') 


dr4:TT£H{nB)r^ 


dip 

dr 


c2 \dr'J 


+ U{R)\ . 


(96) 


Here U{R) is given by Eq. (92), and ip{r) is the velocity potential derived from the wave equation (31) and the 
boundary condition (32) in which ni = rib^H and n 2 = ‘nb,q- 


^In the present analysis, no attention is paid to the dynamical compressibility of the deconfined phase, which may affect the 
kinetic term included in the effective action through a change in the boundary condition at the droplet surface having the form 
(32). 
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B. Nucleation time 


Let us now proceed to evaluate the time r needed to form a real quark matter droplet via quantum tunneling 
in the metastable phase of /3-equilibrated hadronic matter with (without) hyperons, in the regime of overpressures 
and relatively low temperatures, by using the effective action (90) [(96)]. We may obtain t within the exponential 
accuracy as 

r = ly-iexp(^) , (97) 

where t'o is the frequency of the zero-point oscillation around i? = 0, and A is the corresponding extremal value 
of the effective action SeS- Here we have ignored the thermal effect on the pre-exponential factor of r, since the 
frequencies of the higher order oscillations are logarithmically as large as vq. It may also be assumed that the energy 
dissipation and dynamical compressibility in hadronic matter have no influence on the pre-exponential factor, because 
the velocities vf^h and Cg take on nearly the same order of magnitude as the rate |i?| of growth of a virtual droplet. 
The values of obtained by using Eq. (8) amounts to ~ s, being the time scale for strong interactions, 

so we set = 10“^^ s. 

In the absence of hyperons, we recall that the extremal effective action for the quantum nucleation may be roughly 
expressed as Eq. (91). Since the quantity s{T) included in Eq. (29) depends only weakly on temperatures up to 
the quantum-thermal crossover value that will be found to be ^ 10 MeV, we may estimate the time r required to 
form a single quark matter droplet by setting the temperature to be zero. The zero-temperature values of r obtained 
for b = 100 MeV fm“^, ag = 0.4, nig = 0, cTg = 10,30,50 MeV fm“^, and a = 1 have been depicted in Fig. 7(a) 
as a function of P. With the increment in P from the pressure, Pq = 453 MeV fm“^, of the static deconfinement 
transition, the critical droplet radius Rc decreases as Fig. 6(a) exhibits. Such behavior leads to a decrease in A (oc R^) 
and hence to an exponential reduction of t. This indicates that at fixed P, t decreases exponentially as crs(oc Rc) 
decreases linearly. As a consequence, the ambiguity in ag disperses the overpressure AP = P — Pq between ~ 30 and 
~ 140 MeV fm“^ at a typical r. For comparison, the formation time r in the dissipationless case (a = 0) has also 
been evaluated by using the formalism written in Sec. II A; the results have been plotted in Fig. 7(a). The role of 
energy dissipation in increasing the overpressure AP needed to give a constant r has been thus clarified. 

In the case in which the hadronic phase contains hyperons, we may roughly estimate the extremal effective action 
by using Eq. (96). Recall that the extremal value of the compressibility term Sc is likely to behave as oc (T/Tq)^ 
for T <C To, as stated in Sec. II D. Here To is the quantum-thermal crossover temperature found from Eq. (26) as 
A{Eq) = Po/To, where A{Eo) is obtained from Eq. (18) at a given P and T = 0. Since the value of Tq that will be 
estimated to be > 10 MeV is far larger than the temperature of matter in the neutron star core Tns ^0.1 MeV, 
we ignore the contribution of the dynamical compressibility of the hadronic phase to the extremal effective action. 
Thermal effects excite a droplet oscillating around P = 0 from the ground state of energy Eq to the mth bound 
state of energy Em according to the Boltzmann distribution, and hence produce a probability of quantum-mechanical 
formation of a real droplet having the energy Em, relative to the ground-state tunneling probability, as 

Pm = exp j, exp[A(Po) - A{Em)] ■ (98) 

Here to = 1, 2, 3,..., and Em is determined from I{Em) = 2'!T{m+mo + i/A), where I{E) and toq are given by Eqs. (16) 
and (17), respectively. At typical pressures, the difference Em — Em-i is of order 100 MeV, whereas A{Em-i) — A{Em) 
takes on a value of order 10. At T = TnSj the probability Pm is virtually zero, so we take no account of the thermal 
enhancement of the tunneling probability. 

The nucleation time t has been then calculated at zero temperature for hadronic matter with hyperons by substi¬ 
tuting the potential (92) and the effective mass (95) into Eq. (18). In Fig. 7(b) we have plotted the results evaluated 
for b = 100 MeV fm“^, ag = 0.4, nig = 0, and ag = 10,30,50 MeV fm“^, together with the nonrelativistic results 
obtained from Eq. (10). With increasing P from the pressure, Pq = 163 MeV fm“^, of the static deconfinement tran¬ 
sition, the potential barrier is lowered and narrowed for fixed ag, as illustrated in Fig. 6(b). The time r thus shows 
an exponential P dependence in the overpressure regime, as Fig. 7(b) displays. We see by comparison that for large 
ag, the effect of relativity manifests itself as an exponential reduction of t at fixed P. This feature reflects the fact 
that the lower bound Pmin of the positive energy region, given by the maximum value of the quantity U{R) — 2M{R) 
behaving as Fig. 6(b) exhibits, renders the ratio Eq/Uq higher than its nonrelativistic value. The larger ag, the more 
relativistically the droplet moves under the potential barrier. Accordingly, the extent of the reduction in r develops 
at constant P, as can be seen from Fig. 7(b). For small ag, on the other hand, the nonrelativistic results for r are 
slightly smaller than the relativistic results. This is because relativistic corrections act to enhance the kinetic energy 
of the droplet, as is evident from Eq. (12). We also observe in Fig. 7(b) that due to the uncertainty in ag, AP spreads 
between ~ 10 and ~ 60 MeV fm“^ at a typical t. 
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C. stellar conditions 


We turn to the evaluation of the time needed to form the first quark matter droplet in a neutron star whose core 
consists of /3-stable hadronic matter with and without hyperons. This evaluation may be made according to t/N^, 
where Ng is the number of virtual centers of droplet formation in the star, and r is the formation time in the hadronic 
system including only a single nucleation center as obtained in the preceding subsection. Here we have assumed the 
nucleation to be homogeneous: The star retains its hydrostatic equilibrium configuration during its evolution and 
contains no impurities leading to the formation of a seed of quark matter. Ng can be estimated as 


Ng 


3Vc 

4:TTR^ 


(99) 


where 14 is the volume of the central region in which the value of r determined by its pressure remains within an 
order of magnitude of the value calculated at the star’s center. The volume 14 has been obtained as 14 ~ 10® m®, 
irrespective of the presence of hyperons, from the results for the pressure and density profiles of the star that were 
calculated in Sec. Ill from the Tolman-Oppenheimer-Volkoff equation and the equation of state of hadronic matter. 
Hereafter, we set Ng = 10^® since Rc ranges typically ~2-15 fm. The uncertainty in Ng, expected to be an order of 
magnitude or so, has little consequence. 

In order that the first droplet may form during a time scale for the star’s spin-down or accretion Tg ~ yr-Gyr, the 
overpressure should amount to such a value APc as to make r/Ng comparable to Tg. For example, we consider the 
case in which b = 100 MeV fm“®, Ug = 0.4, mg = 0, and a = 1 as exemplified in Fig. 7. For hadronic matter with 
(without) hyperons, the critical overpressure APc lies between ~10 (~30) and ~60 (~140) MeV fm“® because of the 
uncertainty in the surface tension. Correspondingly, the critical droplet radius Rc and the potential barrier height Uq 
take on a value of ~9-10 (~2.0-2.2) fm and ~2-8xl0® (~100-400) MeV, respectively. 

By using the barrier height C/q thus evaluated, we have calculated the quantum-thermal crossover temperature 
To from the relation A = Uq/Tq-, the results have been plotted in Fig. 8 . For comparison, we have also depicted 
the results obtained in the cases of strong dissipation a = 0.2, 5 and no dissipation a = 0 for nuclear matter as 
well as in the nonrelativistic case for hyperonic matter. We see from Fig. 8 (a) that in the absence of hyperons, the 
crossover temperature Tq is large compared with the temperature Tns of matter in the stellar core, although the 
dissipation effects decrease Tq considerably. This result suggests that the formation of a quark matter droplet in 
the core composed of nuclear matter, if occurring, is likely to proceed from quantum fluctuations rather than from 
thermal fluctuations, as is consistent with the result of Ref. [12]. Figure 8 (b) shows that in the presence of hyperons, 
the obtained values of Tq are again large relative to Tns, and that the effect of relativity helps to enhance Tq for large 
ag. To be noted, however, is that the present estimate of Tq (^ 10 MeV) is crude in itself. For such temperatures, 
the thermal effects may alter the thermodynamic quantities included in the effective action (96), induce the heat 
conduction by frequent collisions between the excitations, enhance the rate of quantum-mechanical formation of a 
droplet by raising its energy level around R = 0 from the ground state to the low-lying excited states, and render the 
dynamical compressibility of the hadronic phase effective at increasing the formation rate. At T = Tns, these effects 
are of little consequence, so it is safe to conclude Tq Tns- We have also confirmed that in the present calculations of 
To allowing for and ignoring the presence of hyperons, the order of magnitude of Tq is invariant over the adopted range 
of the bag-model parameters. We remark in passing that a difference in r/Ng by 10 orders of magnitude, expected 
from the ambiguities in Tg, Ng, and vq, alters APc, Rc, and Tq by less than 10%. 

We finally search for the range of the bag-model parameters b, ag, and mg in which the quantum-tunneling nucleation 
of quark matter is expected to occur in the star during the compression of hadronic matter due to the stellar spin-down 
or accretion. Let us now consider the critical condition of the bag-model parameters for the presence of quark matter 
in the star: Initially, the star containing a hadronic core has a central pressure lower than Tmax, and subsequent 
compression, which increases the central pressure up to Pmax in a time scale of Tg, drives the nucleation of quark 
matter in the central region. Here we do not inquire whether the resulting star continues to be gravitationally stable 
or not. The combinations of b and Ug have been then determined for hadronic matter with and without hyperons in 
such a way as to satisfy the critical condition. In the absence of hyperons, we take a = 1 and ag = 30 MeV fm“^; in 
the presence of hyperons, mg and ag have been set as mg = 0,150, 300 MeV and ag = 30 MeV fm“^. The results have 
been plotted in Fig. 4. In these four cases, the critical overpressure APc required to form the first droplet in the star 
results in the shrinkage of the parameter range appropriate to the presence of quark matter from that surrounded by 
the contour of Pq = Tmax- We observe in the present calculations allowing for the presence of hyperons that a peak 
structure appears in the boundary of such a parameter range. Along this boundary, the baryon density discontinuity 
'nb,Q — nb,H at deconfinement runs from ~ —0.1 to ~ 0.3 fm“® with changing b from 50 to 200 MeV fm“®. On the way, 
therefore, the point at which nb,Q = nb,H does exist. At this point, where the effective droplet mass M{R) given by 
Eq. (95) vanishes and so does the critical overpressure APc, the boundary touches the contour of Tq = T^ax, leading 
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to the peak structure. It is noteworthy that around the peak where the effective mass M{Rc) of the critical droplet 
is comparable to or smaller than the barrier height Uq, the effect of relativity on the quantum nucleation becomes 
prominent. 


V. CONCLUSIONS 

We have examined the question of how the presence of hyperons in cold dense hadronic matter affects the dynamical 
properties of the first-order deconfinement transition that may occur in neutron stars. For this purpose, we have used 
the bag-model description of quark matter and the relativistic-mean-field description of hadronic matter. The time 
needed to form a quark matter droplet has been then calculated by using the quantum-tunneling nucleation theory 
which allows for the electrostatic energy and the effects of relativity, energy dissipation, and finite compressibility. 

In the calculations ignoring the presence of hyperons, we have found that the effective mass of a virtual droplet 
moving under the potential barrier is sufficiently large to keep its growth rate fairly small compared with the Fermi 
velocity averaged over nucleon species. This results in an eminent property of the quantum-tunneling formation of 
a real droplet: The quantum tunneling is controlled by the dissipative processes, in which the nucleon quasiparticles 
moving ballistically collide with the droplet surface, rather than by the reversible underbarrier motion of a droplet 
being in a coherent state. We have likewise found that inside a droplet of relatively small radius, an appreciable 
electric charge arises from the pressure equilibrium between quarks and nucleons immersed in a roughly uniform sea 
of electrons and muons. With increase in the droplet radius, the degenerate gas of electrons and muons becomes 
effective at screening the droplet charge, and at last prevents the Coulomb potential barrier from forming. This 
result indicates that a real droplet, if appearing, would develop into bulk matter. We note that these conclusions 
are consistent with the results of Refs. [11] and [12] obtained using a specific set of the bag-model parameters and a 
rather simplified model for nuclear matter. 

The estimates allowing for the presence of hyperons have shown that the density discontinuity at deconfinement 
and the lepton fraction in the hadronic phase fall in magnitude considerably from the values obtained ignoring the 
presence of hyperons. The resulting lepton fraction is small enough for the droplet charge to be nearly zero, allowing 
us to assume the charge neutrality everywhere in the system. The decrease in the effective droplet mass due to the 
lessened density discontinuity, when sufficient, drives the droplet to move relativistically under the potential barrier. 
Not only does such a droplet motion cause the decelerating effect of energy dissipation on the droplet formation to 
decline, but also necessitates the relativistic treatment of the quantum nucleation theory. A natural extension to the 
relativistic regime has been made, leading to an exponential increment in the tunneling probability. 

Irrespective of whether hyperons exist or not, the finite compressibility of the hadronic phase is unlikely to have 
consequence in decreasing the kinetic energy of a droplet via sound emission at low temperatures relevant to neutron 
star cores. Such temperatures have been found to be small compared with the temperature at which an abrupt 
transition from the quantum-tunneling to the thermal-activation nucleation takes place at a typical overpressure. 

We have lastly investigated the range of the parameters b and as where quark matter may appear in the core of 
a neutron star during the compression of hadronic matter due to the stellar spin-down or accretion. It has been 
shown in the analysis of the static deconfinement transition that for the adopted models for bulk hadronic matter 
and quark matter, the existence of hyperons and the finiteness of iris 200 MeV) do not greatly alter the favorable 
parameter range for the appearance of quark matter in the star. We have also found from the dynamical properties of 
the deconfinement transition that the critical overpressure APc required to form the first droplet in the star narrows 
the favorable parameter range obtained from the static properties into the lower b and lower as regime. The critical 
overpressure APc, however, is quite uncertain, because it depends on the poorly known quark-hadron interfacial 
properties and the unclear interactions of quasiparticles in the medium with the droplet surface. In order to make 
a further investigation, we should take account of the possible effect of nucleon superfluidity, which acts to decrease 
the number of available quasiparticle states with momenta close to the Fermi surfaces [12]. If the nucleons were in 
a superfluid state, the friction force exerted by the quasiparticles on the droplet and hence the critical overpressure 
APc in the system not containing hyperons would be reduced considerably. 

It is important to bear in mind that the obtained results for the deconfinement transition are based on the adopted 
models for hadronic matter and for quark matter, which are inequivalent physical descriptions of the two phases and 
hence do not necessarily describe the transition well. Granted that such models hold good, the dynamical properties 
of the transition are not obvious quantitatively as mentioned in the preceding paragraph. Moreover, the mechanical 
or chemical impurities leading to the formation of a quark matter seed, if occurring in the star at all, might relax the 
condition for the transition drastically. 

Once the first droplet of quark matter arises from fluctuations in a neutron star, its growth into bulk matter would 
cause strange matter to occupy the central region of the star as long as the strange matter is more favorable than 
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/3-equilibrated hadronic matter. These two phases, if both are present, would have a charged but overall neutral layer 
which spreads from the quark-hadron boundary on each side over a length scale determined by the Thomas-Fermi 
screening length. At the boundary, they not only would have the same chemical potential of leptons, but also would 
be in equilibrium with respect to (3 and deconfinement processes. Even if the quark-hadron mixed phase is the lowest 
energy configuration at pressures around the boundary, its appearance depends on whether the boundary layer is 
stable against fragmentation such as the clustering of protons and of s quarks. The elucidation of this problem is 
beyond the scope of this paper. 

We believe that the present analysis is applicable to the study of the kinetics of the deconfinement transition that 
may occur in hot dense matter encountered in stellar collapse and in relativistic heavy-ion collisions. Due to its high 
temperatures (above ~ 10 MeV), the nucleation of quark matter in hadronic matter may proceed classically, i.e., via 
thermal activation. Otherwise, thermally assisted quantum nucleation may take place; it may be influenced by the 
irreversible transport of heat and momentum in the medium and by the dynamical compressibility of the hadronic 
phase. In the case of stellar collapse, neutrinos trapped in a supernova core, leading to large fractions of electrons 
and muons, may cause a quark matter droplet to have an appreciable electric charge. 
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APPENDIX: DISTRIBUTION INSIDE A QUARK MATTER DROPLET 


A quark matter droplet, if having a macroscopic baryon number and a nonzero electric charge in itself, has the 
constituent quarks distributed over a scale of the screening length from the surface just like a conductor. In the case 
of a droplet of /3-stable quark matter in vacuum, such a screening action by quarks was investigated by Heiselberg [40] 
within the linear Thomas-Fermi approximation [35]. We now extend his calculations to the case of a quark matter 
droplet formed in /3-equilibrated nuclear matter under the flavor conservation, considered in Sec. IV A 1, and ask to 
what extent the resulting quark distributions alter the potential for the formation and growth of the droplet. 

We begin with the potential energy U{R) of a droplet of radius R given by Eq. (68), in which we remove the 
assumption that the quark number densities Uq {q = u,d) are constant inside the droplet, and we set 

{ nqo + 5nq{r) , r < R , 

(Al) 

0 , r > R . 

Here Uuo and Udo, independent of r, satisfy the pressure equilibrium condition (76) and the flavor conservation (61). 
In determining the deviations Suq on a footing equal to 6ni, we first note that the baryon density and the electric 
charge density po for the inhomogeneous phase are written as 

nb,D{r) = ^[n„(r) -I- nd{r)] -|- 0{r - R)nbM (A2) 

and 


Pair) = e 


qin^{r) + 9{r - R)np 

i=q,l 


= 9{R - r)pQo + Sp{r) . 


Here J/j is the deviation of the charge density from that calculated in the incompressible limit 


(A3) 


Sp(r) = e 


dq^nq + qi5ni 

Q I 


(A4) 


and pqq is given by Eq. (74) in which n„ = n„ 0 : nd = Udo, and Ug = 0. We then expand the potential U{R) up to 
second order in Suq and Sni , and minimize the resulting potential with respect to Suq and 6ni under the global flavor 
conservation 
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[ dVifoSuu - Sud) = 0 , 

Jr<R 


with /o = the global neutrality of electric charge 


/ dVqiSni + / dVqqSuq + Zq = 0 , 

, Jv ^ Jr<R 


I 

where Zq is given by Eq. (73). This minimization leads not only to Eq. (78) for I leptons, but also to 

^2 

_ I ^ _L // ^ _ j1i »y, _L ^ ,, _1_ / 


K 


47r(gqe)2 


' qqe(j){r) + ^qO - + qq^J^e,H + Cq 


for q quarks, where /i^o is the chemical potential of q quarks calculated at density Uqo, 

1/2 


Kn = 


Is 


Xq — XlqQ 


(A5) 


(A6) 


(A7) 


(A8) 


is the reciprocal of the Thomas-Fermi screening length of q quarks, and C„ = Cfo and Cd = —C are the constants 
calculated at a given R from the flavor conservation (A5). Here the proton screening continues to be ignored. We 
also neglect the energy increment yielded by the gradient of Suq present over a scale of the screening length from 
the droplet surface. This increment, acting to flatten the quark distributions, is essentially unknown along with the 
quark-hadron interfacial property. 

The electrostatic potential obeys the Poisson equation 


V^(/)(r) — n{rY(j>{r) = —47rpQ0(i? — r) , 


where 


PQ = PQO-Z' 


Airq e ^Pb,H + qqPe,H + Cq 


(A9) 


(AlO) 


and 


.) = 


K< = ^/kI + kI + kI , r < R , 
kl , r > R . 


(All) 


Here k< (kl) corresponds to the reciprocal of the screening length inside (outside) the droplet. The solution to Eq. 
(A9) is obtained as 


47rpQ 


(/)(r) = < 


1 - 


K<i?(l -I- KlR) 


sinh(K<r) 


KLi?sinh(K<i?) -I- K<i?COsh(K<i?) K<r 


, r < R 


KLi?[K<i?cosh(/t<i?) — sinh(K<i?)] exp[Ki(i? — r)] 
[ K% Ki,7?sinh(K<i?)-I-K<i?cosh(K<i?) klt 


r > R . 


It is straightforward to derive the electrostatic energy Ec from Eq. (A12); the result is 

^2 J2 


EciR) = ^ 


1 


[KLi?sinh(K<i?)-I-K<i?COsh(K<i?)]^ 

2sinh^(K<i?) 


K<i? -I- cosh(K<i?) sinh(K<i?) — 


K<i? 


{K<Rf{l + KLRf 


K< 

kl 


[K<i?COSh(K<i?) — sinh{K^R)]^KLR{KLR + 2) 


(A12) 


(A13) 


The quark distributions replace the excess lepton energy AEl{R) in the potential (68) with the sum of AEl{R) and 
the excess quark energy defined as 
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(A14) 


AEq{R) = / dV{eQ[nq{r)] - sgiugo)} . 

Jr<R 

The linear approximation used here yields 

AEl{R) + AEq{R) = ZoHe,H + y] ^^90 - J SUgEV . (A15) 

By using Eqs. (78), (A4), and (A7), we have calculated the density deviations Suq, 6ni, and Sp/e from the values 
obtained not allowing for the quark or lepton screening. We have confirmed that for i? < 10 fm, the number density 
of each component remains close to its constant value obtained in the incompressible limit. The results calculated for 
the inhomogeneous phase characterized by 6 = 100 MeV fm“^, as = 0.4, i? = 2,5 fm, and P = 525 MeV fm“^ have 
been plotted in Fig. 9, together with the charge density deviation calculated allowing for the lepton screening alone 
as in the main text. We see from the comparison of these charge density deviations Sp that the quark screening does 
not vary the total charge distribution pD significantly even for i? = 5 fm, although ~ 3 fm is small relative to this 
R and to ~ 5 fm. This result reflects the feature that under the flavor conservation, the static compressibility 
of the quark fluid plays a role in chemically decreasing the sum of the excess lepton and quark energies rather than 
in reducing the electrostatic energy via Coulomb screening. Note that 6nq > 0, leading to a reduction of the sum 
AEl + AEq given by Eq. (A15). The resulting decrease in the potential U{R) can be observed in Fig. 10 by comparing 
the result obtained for a lepton- and quark-screened droplet with that for a lepton-screened droplet. The case of a 
lepton- and quark-screened droplet agrees fairly well with the case of a lepton-screened droplet and the incompressible 
case in which we set Eq = 0, since the inhomogeneous phase generally satisfies 

ZoPe,H > Zop,e,H “ {AEl + AEq) > Eq > \AEc\ , 

where Ec is given by Eq. (A13), and AEc is the difference between the right-hand sides of Eqs. (80) and (A13). It 
is of interest to note that the finite compressibility of the fluid of N nucleons produces no chemically induced change 
in their number density, as can be found from the expansion of U{R) up to second order in Sun- This property is 
due to the {3 equilibrium sustained in the hadronic phase. 

The shifts (5n„ and 6nd due to the quark screening, being typically ~ 0.1 fm“^, in turn affect the kinetic properties 
of the droplet. First, the effective droplet mass M{R), estimated from the nucleon kinetic energy as Eq. (85), is 
increased, since it is proportional to the square of the baryon density discontinuity &t r = R. Such increase in 
M{R) diminishes the typical velocity of the precritical droplet, enforcing its smallness relative to the averaged Fermi 
velocity fF.ff |nB^iv=o given by Eq. (89). On the other hand, the hydrodynamic mass flow of each quark component, 
which arises from the Coulomb-induced inhomogeneity in Sug, contributes little to M{R). Second, the friction force 
FH\nB^N=o given by Eq. (87), behaving as EH\nB^N=o oc M{R), likewise becomes strong owing to Sug. A further 
friction force arises from the collisions of quark quasiparticles, carrying the quark mass flow, with the droplet surface, 
but this is negligibly small compared with the nucleon contribution. Recall that the baryon density discontinuity at 
r = R is fairly large, ranging ~ 0.1-0.7 fm“^, in the absence of the quark screening. We may then conclude that the 
shifts Suq, acting to enhance the dissipation effects governed by nucleon quasiparticles, lead to a fractional increase 
in the critical overpressure APc- 

We conclude this appendix by briefly mentioning the influence of the static compressibility of quark matter, sur¬ 
rounded by hadronic matter with hyperons, on the potential energy. This effect can easily be estimated by following a 
line of argument analogous to the case in which no hyperons exist and by neglecting the electrostatic properties. Such 
an estimate made over the considered range of the parameters such as b, as, Us, and rus demonstrates that the finite 
compressibility increases the number densities n„, Ud, and n*, independently of R, by less than about 10% at typical 
overpressures. This leads to a chemically induced decrease in the quark energy, which has been found to be of little 
consequence to the potential energy. The resulting increase in the baryon density inside a droplet by typically ~ 0.05 
fm“^ shifts the baryon kinetic energy appreciably, and hence the effective droplet mass (95). This is because, in the 
absence of the quark screening, the baryon density discontinuity at r = i? is relatively small, ranging ~ —0.1-0.3 
fm“^. Such shifts move into the lower-6 regime the region of the bag-model parameters in which the droplet be¬ 
haves relativistically, but otherwise do not change the qualitative conclusions about the nucleation processes obtained 
ignoring the quark screening. 
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FIG. 2. Composition of /J-stable nuclear matter (upper panel) and hyperonic matter (lower panel) at zero temperature, 
calculated using the model of Glendenning and Moszkowski [16]. 
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FIG. 3. Equation of state of zero-temperature nuclear matter in f) equilibrium. The solid and long-dashed lines are the 
results calculated from the model of Glendenning and Moszkowski [16] within the relativistic mean-field (RMF) theory; the 
latter has been obtained allowing for the presence of hyperons. The dash-dotted and dotted lines are the results obtained by 
Wiringa, Fiks, and Fabrocini [25] using the Urbana Vi 4 (UV14) two-nucleon potential plus the Urbana VII (UVII) three-nucleon 
potential and the density-dependent three-nucleon interaction (TNI) model, respectively. The short-dashed line is the relativistic 
Brueckner-Hartree-Fock (RBHF) result obtained by Engvik et al. [26] using the Bonn A two-nucleon potential. 
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FIG. 4. Contour plots of the pressure Po (as marked in MeV fm“®) of the static deconfinement transition on the b-Os plane. 
The solid lines in (a) denote the contours of the transition pressures from /3-stable nuclear matter to u and d quark matter. The 
solid lines in (b), (c), and (d) represent the contours of the transition pressures from /3-stable hyperonic matter to u, d, and s 
quark matter calculated for ms = 0,150, 300 MeV, respectively. In each panel, the contours of Pq = Pmax (short-dashed line), 
Po = Pi. 4 (dotted line), and Po = P^- (dash-dotted line) are included (see text). We have also plotted the combinations of b 
and as satisfying the critical condition required to form a real quark matter droplet via deconhnement in the star that contains 
the nucleation centers of number Ns = 10^® and is spinning down or accreting matter from a neighboring star in a time scale 
of Ts = 1 Myr. The lower b and lower a® region is favorable for the presence of quark matter in the star. The long-dashed 
line in (a) denotes the result obtained for a droplet arising from nuclear matter by setting as = 30 MeV fm“^ and a = 1; the 
long-dashed line in (b)--(d), for a droplet deconflned from hyperonic matter by setting as = 30 MeV fm“^. 
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FIG. 5. (a) [(b)] Chemical potentials for the electrically neutral bulk phases as a function of pressure, calculated for the 
bag-model parameters b — 100 MeV fm“®, Os = 0.4, and ms = 0. The solid line represents the phase of /3-equilibrated hadronic 
matter without (with) hyperons; the dashed line, the phase of matter deconfined therefrom; the dotted line, the phase of strange 
matter, i.e., u, d, and s quark matter in /3 equilibrium, (c) [(d)] Energy densities for the corresponding phases as a function 
of baryon density. The dash-dotted lines are the double-tangent constructions denoting the coexistence of the two bulk phases 
involved. 
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FIG. 6. Potential energies U{R) of a quark matter droplet present in /3-stable nuclear matter (a) and hyperonic matter (b), 
evaluated for the bag-model parameters b = 100 MeV fm“®, as = 0.4, and = 0 and for various pressures above the pressure 
Po of the static deconfinement transition. In (a), we have set as = 30 MeV fm“^. The solid lines are the results obtained 
from Eq. (68) for a lepton-screened droplet; the dashed lines, for a nonscreened droplet; the dash-dotted lines, in the case in 
which we set Ec = 0. The solid lines in (b) are the results evaluated from Eq. (92) for as — 30 MeV fm“^. The quantities 
U{R) — 2M{R), where M{R) is the effective droplet mass given by Eq. (95), have also been plotted in (b) by dashed lines. The 
maximum value of U(R) — 2M{R) yields the upper bound of the forbidden region of the droplet energy. 


31 









150 200 250 300 

P (MeV fm-3) 


FIG. 7. Time needed to form a real quark matter droplet via quantum tunneling in /I-stable zero-temperature nuclear matter 
(a) and hyperonic matter (b), evaluated as a function of pressure for the bag-model parameters b — 100 MeV fm“®, as = 0.4, 
and nis = 0. In (a), the surface tension has been set as Gs = 10, 30, 50 MeV fm“^. The solid lines are the dissipative results 
obtained for a = 1, and the dashed lines are the nondissipative results. In (b), the relativistic (solid lines) and nonrelativistic 
(dashed lines) results calculated for as = 10, 30, 50 MeV fm“^ are plotted. For comparison, we have marked in both panels the 
pressure Po of the static deconfinement transition and the central pressure Pmax of the maximum-mass neutron star. 
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FIG. 8. Crossover temperatures To from the quantum-tunneling to the thermal-activation nucleation of quark matter in 
/3-stable nuclear matter (a) and hyperonic matter (b), calculated as a function of the surface tension for the bag-model 
parameters b — 100 MeV fm“®, as = 0.4, and rus = 0 and for the stellar conditions t/Ns = 1 Myr and Ns = 10^®. In (a), the 
solid lines are the results obtained for dissipative cases of a = 0.2,1, 5, and the dashed line is the nondissipative result. In (b), 
the solid and dashed lines are the relativistic and nonrelativistic results, respectively, estimated not allowing for thermal effects 
on the quantum nucleation. 


33 






0 2 4 6 8 10 


r (fm) 

FIG. 9. Deviations 5n„, Sud, and drie+^n^i of the quark and lepton number densities from the values n-uo, ndo, and ne,ir+n^,_tf 
obtained in the incompressible limit due to the screening (static compressibility) of the fluids of quarks and leptons, calculated 
ignoring the presence of hyperons for the droplet radius R — 2 fm (upper panel) and i? = 5 fm (lower panel), the pressure 
P = 525 MeV fm“®, and the bag-model parameters b = 100 MeV fm“® and a* = 0.4. The resulting deviation 5p of the electric 
charge density from the stepwise charge distribution obtained in the incompressible limit has also been plotted, together with 
the result allowing for the lepton screening alone (no quark screening). Since Sue graphically agrees with only the sum 
Srie + Sri/j, has been shown. See the text for the definition of the screening lengths and 
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FIG. 10. Potential energies U{R) of a qnark matter droplet in the metastable phase of /3-stable nuclear matter evaluated 
for the surface tension as = 30 MeV fm~^, the bag-model parameters b = 100 MeV fm“® and as = 0.4, and various pressures 
above the pressure, Po « 453 MeV fm“®, of the static deconfinement transition. The solid lines are the results obtained from 
Eq. (68) for a lepton- and quark-screened droplet; the long-dashed lines, for a lepton-screened droplet; the short-dashed lines, 
for a nonscreened droplet; the dash-dotted lines, in the incompressible case in which we set Ec = 0. 
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TABLE I. Values of M^ax, ^max, -Pi. 4 , and P^- evaluated for hyperonic and nuclear matter using the relativistic mean-field 
theory with the parameters taken from Ref. [16]. 




P 

max 

P1.4 

P^- 


A^max 

(MeV fm“®) 

(MeV fm“®) 

(MeV fm“®) 

Hyperonic matter 

1.78M0 

261.0 

43.13 

23.30 

Nuclear matter 

2.36M0 

500.4 

38.91 
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